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1. Introduction 

In this paper, we study the convergence of certain nonlinear iterative reconstruction 
methods for inverse problems in Banach spaces. We consider a class of inverse problems 
defined by a nonlinear map from parameter or model functions to the data. The 
parameter functions and data are contained in certain Banach spaces, or Hilbert 
spaces, respectively. We explicitly construct sequences of parameter functions by a 
Landweber iteration. Our analysis pertains to obtaining natural conditions for the 
strong convergence of these sequences (locally) to the solutions in an appropriate 
distance measure. 

Our main result establishes convergence of the Landweber iteration if the inverse 
problem ensures a Holder stability estimate. Moreover, we prove monotonicity of 
the residuals defined by the sequence induced by the iteration. We also obtain the 
convergence rates without so-called source and nonlinearity conditions. The stability 
condition is a natural one in the framework of iterative reconstruction. 

Extensive research has been carried out to study convergence of the Landweber 
iteration [25] and its modifications. In the case of model and data spaces being Hilbert, 
see Hanke, Neubauer & Scherzer [18]. An overview of iterative methods for inverse 
problems in Hilbert spaces can be found, for example, in Kaltenbacher, Neubauer & 
Scherzer [21]. Schopfer, Louis & Schuster [30] presented a nonlinear extension of the 
Landweber method to Banach spaces using duality mappings. We use this iterative 
method in the analysis presented here. Duality mappings also play a role in iterative 
schemes for monotone and accretive operators (see Alber [2], Chidume & Zegeye [H] 
and Zeidler [36| [37]). The model space needs to be smooth and uniformly convex, 
however, the data space can be an arbitrary Banach space. Due to the geometrical 
characteristics of Banach spaces other than Hilbert spaces, it is more appropriate to use 
Bregman distances rather than Ljapunov functionals to prove convergence (Osher et al. 
[27]). For convergence rates, see Hofmann et al. [20]. Schopfer, Louis & Schuster [3~T] 
furthermore considered the solution of convex split feasibility problems in Banach spaces 
by cyclic projections. Under the so-called tangential cone condition, pertaining to the 
nonlinear map modelling the data, convergence has been established; invoking a source 
condition in a convergence rate result. Here, we build on the work of Kaltenbacher, 
Schopfer and Schuster [22] and revisit these conditions with a view to stability properties 
of the inverse problem. 

In many inverse problems one probes a medium, or an obstacle, with a particular 
type of field and measures the response. From these measurements one aims to 
determine the medium properties and/or (geometrical) structure. Typically, the physical 
phenomenon is modeled by partial differential equations and the medium properties 
by variable, and possibly singular, coefficients. The interaction of fields is usually 
restricted to a bounded domain with boundary. Experiments can be carried out on 
the boundary. The goal is thus to infer information on the coefficients in the interior of 
the domain from the associated boundary measurements. The map, solving the partial 
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differential equations, from coefficients or parameter functions to the measurements or 
data is nonlinear. Its injectivity is studied in the analysis of inverse problems. As an 
example, we discuss Electrical Impedance Tomography, where the Dirichlet-to-Neumann 
map represents the data, and summarize the conditions leading to Lipschitz stability. 

Traditionally, the Landweber iteration has been viewed as a fixed-point iteration. 
However, in general, for inverse problems, the underlying fixed point operator is not 
a contraction. There is an extensive literature of iterative methods for approximating 
fixed points of non-expansive operators. Hanke, Neubauer & Scherzer [18] replace the 
condition of non-expansive to a local tangential cone condition, which guarantees a local 
result. In the finite-dimensional setting, in which, for example the model space is M. n , 
non-convex constraint optimization problems admitting iterative solutions have been 
studied by Curtis et al. [TB] . Under certain assumptions, they obtain convergence to 
stationary points of the associated feasibility problem. In the context of inverse problems 
defined by partial differential equations, this setting is motivated by discretizing the 
problems prior to studying the convergence (locally) of the iterations. Inequality 
constraints are necessary to enforce locality. The non-convexity is addressed by Hessian 
modifications based on inertia tests. 

The paper is organized as follows. In the next section, we summarize certain 
geometrical aspects of Banach spaces, including (uniform) smoothness and (uniform) 
convexity, and their connection to duality mappings. Smoothness is naturally related to 
Gateaux differentiability. We also introduce the Bregman distance. We then define the 
nonlinear Landweber iteration in Banach spaces. In Section [3] we introduce the basic 
assumptions including Holder stability and analyze the convergence of the Landweber 
iteration in Hilbert spaces. In Section H] we adapt these assumptions and generalize the 
analysis of convergence of the Landweber iteration to Banach spaces. We also establish 
the convergence rates. In Section [5] we give an example, namely, the reconstruction of 
conductivity in Electrical Impedance Tomography, and show that our assumptions can 
be satisfied. 

2. Landweber iteration in Banach spaces 

Let X and Y be both real Banach spaces. We consider the nonlinear operator equation 



with domain T>(F) C X. In applications, F : T>{F) — > Y models the data. In the 
inverse problem one is concerned with the question whether y determines x. We assume 
that F is continuous, and that F is Frechet differentiate, locally. 

We couple the uniqueness and stability analysis of the inverse problem to a local 
solution construction based on the Landweber iteration. Throughout this paper, we 
assume that the data y in (12.11) is attainable, that is, that (12. ip has a solution x* (which 
needs not to be unique). 





(2.1) 
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2.1. Duality mappings 

The duals of X and Y are denoted by X* and Y* , respectively. Their norms are denoted 
by || • || . We denote the space of continuous linear operators X — > Y by C(X,Y). Let 
A : V(A) C X ->• Y be continuous. Here V(A) denotes the domain of A. Let heV(A) 
and k G X and assume that h + t(/c — /i) G f(^4) for all t G (0, to) for some to, then we 
denote by DA(h)(k) the directional derivative of A at /i G ©(A) in direction k G ^(-4), 
that is, 

DA(h)(k):=^ A(h + t ^- A(h) . 

If DA(h) G £(X, F), then DA{h) is called Gateaux different iable. If, in addition, the 
convergence is uniform for all k G B to , then DA is Frechet differentiate at h. For 
x G X and x* G X*, we write the dual pair as {x,x*} = x*(x). We write A* for the 
dual operator A* G C(Y*,X*) and ||A|| = \\A*\\ for the operator norm of A. We let 
1 < p, q < oo be conjugate exponents, that is, 

- + - = 1. (2.2) 
P 9 

For p > 1, the subdifferential mapping J p = df p : X — > 2 X * of the convex functional 
fp '. x i — y ^||^|| p defined by 

J p (x) = {x* G X* | (x,x*) = \\x\\ \\x*\\ and \\x*\\ = IM^ 1 } (2.3) 

is called the duality mapping of X with gauge function t i->- t p_1 . Generally, the duality 
mapping is set-valued. In order to let J p be single valued, we need to introduce the 
notion of convexity and smoothness of Banach spaces. 
One defines the convexity modulus 5x of X by 

S x (e) — inf {1 - \\Ux + x)\\ ] \\x\\ — \\x\\ — 1 and \\x - x\\ > e} (2.4) 

x,xex 

and the smoothness modulus px of X by 



Pxij) = sup {|(||^ + tx\\ + \\x — rx\\ — 2) | \\x\\ = \\x\\ = 1}. (2.5) 

x,x£X 



Definition 2.1 ^4 Banach space X is said to be 

(a) uniformly convex if 5x{t) > for any e G (0,2], 

(b) uniformly smooth z/lim T _ >0 Px ^ = 0, 

(c) convex of power type p or p-convex if there exists a constant C > such that 
Sx(e) > Ce?, 

(d) smooth of power type q or q-smooth if there exists a constant C > such that 
Px(r) < Ct«. 



Example 2.2 (a) A Hilbert space X is 2-convex and 2-smooth and J<i : X — > X is the 
identity mapping. 
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(b) Let Q C M n be an open domain. The Banach space LP = L P (Q), p > 1 is uniformly 
convex and uniformly smooth, and 

e 2 , l<p<2, 
e p , 2 < p < oo; 



Plv{t) 



t p , 1 < p < 2, 
T 2 , 2 < p < oo. 



(c) For X = L r (R n ), r > 1, we have 
J p : L r (M") -»■ 

u(a;) i — ^ ||n||^7 r |M(x)| r_2 M(a;), 



where - + - = 1. 

For a detailed introduction to the geometry of Banach spaces and the duality 
mapping, we refer to (151 [30] . We list the properties we need here in the following 
theorem. 



Theorem 2.3 Let p > 1. The following statements hold true: 

(a) For every x G X, the set J p (x) is not empty and it is convex and weakly closed in 
X*. 

(b) If a Banach space is uniformly convex, it is reflexive. 

(c) A Banach space X is uniformly convex (resp. uniformly smooth) iff X* is uniformly 
smooth (resp. uniformly convex). 

(d) If a Banach space X is uniformly smooth, J p (x) is single valued for all x G X. 

(e) If a Banach space X is uniformly smooth and uniformly convex, J p (x) is bijective 
and its inverse J" 1 : X* — >■ X is given by J' 1 = J* with J* being the duality 
mapping of X* with gauge function t i-> t 9-1 , where 1 < p, q < oo are conjugate 
exponents. 



Throughout this paper, we assume that X is p-convex and g-smooth with p, q > 1, 
hence it is uniformly smooth and uniformly convex. Furthermore, X is reflexive and 
its dual X* has the same properties. Y is allowed to be an arbitrary Banach space; j p 
will be a single-valued selection of the possibly set-valued duality mapping of Y with 
gauge function t H- t p_1 , p > 1. Further restrictions on X and Y will be indicated in 
the respective theorems below. 
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Figure 1. Brcgman distance, A p . 



2.2. Bregman distances 

Due to the geometrical characteristics of Banach spaces different from those of Hilbert 
spaces, it is often more appropriate to use the Bregman distance instead of the 
conventional norm-based functionals \\x — x\\ p or ||</ p (:e) — J p (x)|| p for convergence 
analysis. This idea goes back to Bregman [TT] . 

Definition 2.4 Let X be a uniformly smooth Banach space and p > 1. The Bregman 
distance A p (x, •) of the convex functional x i— > ^\\x\\ p at x G X is defined as 

A p (x,x) = -\\x\\ p \\x\\ p — (J„(x) , x — x) , x6l, (2.6) 

p p 

where J p denotes the duality mapping of X with gauge function t \-t t p ~ x . 

In the following theorem, we summarize some facts concerning the Bregman 
distance and the relationship between Bregman distance and the norm [Tl l2l [121 [35] . 

Theorem 2.5 Let X be a uniformly smooth and uniformly convex Banach space. Then, 
for all x, x G X , the following holds: 



(a) 



AJx,x) = -\\x\\ p - -\\x\\ p - (JJx),x) + \\x\\ p (2.7) 
p p 

— 1 1 X 1 1 I 1 1 £C 1 1 \ J p (y X ^ . CC ^ • 

p q 



(b) A p (x, x) > and A p (x, x) = -v^ x = x. 

(c) A p is continuous in both arguments. 

(d) The following statements are equivalent 

(i) lim^oo \\x n -x\\ = 0, 
(ii) lim^oo A p (x n ,x) = 0, 
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(Hi) lim n _^ 00 ||x n || — ||*^|| and lim^^oo^i/p^n), x) — (t/p(x),x). 

(e) If X is p-convex, there exists a constant C p > such that 

A„(x,x) > ^||x-xf. (2.8) 
V 

(f) If X* is q-smooth, there exists a constant G q > such that 

AJx*,x*) < -l\\x*-x*\\ q , (2.9) 
Q 

for all x*,x* G X*. 



Remark 2.6 The Bregman distance A p is similar to a metric, but, in general, does not 
satisfy the triangle inequality nor symmetry. In a Hilbert space, A 2 ( 



2.3. Landweber iteration 

In this subsection, we introduce an iterative method for minimizing the functional 

= -\\F(x) -y\\ p . (2.10) 
p 

The iterates {x^} are generated with the steepest descent flow given by 

d& k \x k ) = DF(x k )*j p (F(x k ) - y). (2.11) 
To be more precisely, we study the iterative method in Banach spaces, 
J p (x k+1 ) = J p {x k ) - fj l DF(x k )*j p (F(x k ) - y), 

X k+1 = J*(J p (xk+l)), 

where J p : X — > X*, J* : X* — > X and j p : Y — > Y* denote duality mappings in 
corresponding spaces. When X and Y are Hilbert spaces and p = 2, this reduces to the 
Landweber iteration in Hilbert spaces 

x k+1 = x k - /iDF(x k )*(F(x k ) - y). (2.13) 

If F is a linear operator, the iteration (1.3) coincides with Landweber 's original 
algorithm. We specify \x below. Equation (12 . 1 21) defines a sequence (x k ). 
If F(x*) = y, the so-called tangential cone condition [22J, 

\\F(x)-F(x)-DF(x)(x-x)\\ < c TC ||F(x)-F(£)|| Vx, x G ^(x f ),(2.14) 

for some < ctc < 1, is crucial to obtain convergence of (x k ) to x^ [THJ EUJ [22] : 
Bf( x ^) = {x £ X \ A p (x,x^) < p} C T>(F). A source condition controls the 
convergence rate. Here, we study convergence and convergence rates in relation to a 
single, alternative condition replacing the tangential cone and source conditions, namely, 
Holder type stability, 

A p (x,x) < C p F \\F(x) -F{x)\\ 1 ¥ lp Vx, x G Bf{x ] ), 

for some e G (0, 1]. With the Frechet differentiability of F and the Lipschitz continuity 
of DF, this condition implies the tangential cone condition, and, hence, convergence is 
guaranteed; however, it also implies a certain convergence rate. 
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3. Convergence rate and radius of convergence — Hilbert spaces 

In this section, we assume that X and Y are Hilbert spaces. Then the mappings J p ,j p 
and J* are identity mappings. Let B p (x ) denote a closed ball centered at x with radius 
p, such that B = £y (xo) C T>(F), p' > p. As before, let x* generate the data y, that is 

F(xt)=y. (3.1) 

We assume that x^ G B p (x Q ). 

Assumption 3.1 (a) The Frechet derivative, DF, of F is Lipschitz continuous locally 
in B and 

\\DF(x)-DF(x)\\<L\\x-x\\ Wx,xeB. (3.2) 

(b) F is weakly sequentially closed, that is, 

X n ^X, ]f XEV(F), 

F(x n ) ->y J I F(x) = y. 

(c) The inversion has the uniform Holder type stability, that is, there exists a constant, 
Cf > 0, such that 

-^=\\x-x\\ < C F \\F{x) -Fix)^ \/x,xeB (3.3) 
for some e G (0, 1] 

In the remainder of this section, we discuss the convergence criterion and 
convergence rate for the Landweber iteration (12.131) . 



Theorem 3.2 Assume there exists a solution x^ to 113. 1\) and that Assumption \3. 1\ holds. 
Furthermore, assume that 

\\DF{x)\\<L MxeB. (3.4) 

Let the positive stepsize p be such that 

1 

a < — , 

L\ ^ (3-5) 
p{l - pi 2 ) < 2—C F +s 



Let 



If 



-(2LL £ C 2 F )- 2/£ . 



P 2 



-||ac - 1| 2 < P, (3.6) 
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then the iterates satisfy 



2l|zk-a; + || 2 <P, k — 1,2,... (3.7) 
and it -)• as k —y oo. Moreover, let 

c=^(l-fiL 2 )C~^; (3.8) 
from ( fff. 5\) . it follows that < c < 1 . The convergence rate is given by 

^\\x k -x*\\ 2 <p{l-c) k , (3.9) 
if e = 1. For e G (0, 1), the convergence rate is given by 

-\\x k -xm 2 < (ckj^ + p-^j "\ k = 0,1,.... (3.10) 

The proof is a special case of the Banach space setting, cf. Theorem I4.5( see 
Section HJ The convergence is sublinear if < e < 1 and the speed up as e — > 1 
relates to the fact that it switches to a linear convergence. 

For the critical index e — 0, that is, the power in the right-hand side of the stability 
inequality (13. 3 j) equals to |, we need to invoke an assumption on the stability constant Cp 
to arrive at the convergence and convergence rate results. An interesting by-product is 
that the convergence radius only depends on the radius within which the Holder stability 
(13. 3p holds. Hence, if the forward operator F satisfies (13. 3 j) globally, then we obtain a 
global convergence and convergence rate result. 

Theorem 3.3 Assume there exists a solution x^ to A3. 1\) and that Assumption \3. 1\ holds 
with e = 0. Furthermore, assume that 

\\DF{x)\\<L VxGi3. (3.11) 
Let the stability constant Cp and the positive stepsize fi satisfy that 

pL 2 + 2LC 2 F <2. (3.12) 
Then the iterates satisfy 

Xk — > x^ as k —7- oo. 

Moreover, let 

c = ^(-2 + pL 2 + 2LC 2 p)C- F \ (3.13) 
The convergence rate is given by 

d|z fc - x j \\ 2 < (2\\x - a: 1 "!!" 2 + ck)- 1 . (3.14) 



The proof is again a special case of the Banach space setting, cf. Theorem I4.5j see 
Section HI 
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Remark 3.4 The convergence radius condition ( Iff. 6]) on the starting point Xq may be 
replaced by a convergence radius condition on the starting simulated data F(xq), 

\\F(x )-y\\ 1+£ <pCf. (3-15) 

In fact, with the aid of the stability inequality l\3.3\) , 113.6}) follows from Ii3.15\) . 

4. Convergence rate and radius of convergence — Banach spaces 

In this section, we discuss the convergence and convergence rate of the Landweber 
iteration (12.121) in Banach spaces. Let B p {xq) denote a closed ball centered at xq with 
radius p, and B = B^(x^) denote a ball with respect to the Bregman distance centered 
at some solution x^ . We assume that B^(x^) C T>(F). 

Assumption 4.1 (a) The Frechet derivative, DF , of F is Lipschitz continuous locally 
in B and 

\\DF(x)- DF(x)\\<L\\x-x\\ \/x,xeB. (4.1) 

(b) F is weakly sequentially closed, that is, 

X n ^X, If XE V(F), 

F(x n ) ^y J [ F(x) = y. 

(c) The inversion has the uniform Holder type stability, that is, there exists a constant 
Cp > such that 

A p (x,x) <C p F \\F(x)-F{x)\\ i ¥ p Vx,xeB, (4.2) 

for some e G (0, 1]. 

Remark 4.2 Note that the nonemptyness of the interior (with respect to norm) ofV(F) 
is sufficient for B C T>(F). 



Remark 4.3 With the assumption that X is p-convex, w ^th l[2.8\) implies the 

regular notion of Holder stability in norm. 



Remark 4.4 Under the Lipschitz type stability assumption, that is, lj4-£\ ) with e — I, 
we have that 

(J p (x^), x — x^) < ||:e'I'|| p ~ 1 ||:e — x'\\ 

< CA p (x,it)Vp 

< CC F \\F(x) - F(x^)\\, MxeB 

for some constant C > 0. It has been shown in J2$tf that this implies the source condition, 

J p (x ] ) = DF{x ] )*u 
for some u satisfying < 1. 



Local Analysis of Inverse Problems: Holder Stability and Iterative Reconstruction 11 

Theorem 4.5 Let Y be a general Banach space, and X be a Banach space which is 
p-convex and q-smooth with l/p + 1/q — 1. Assume there exists a solution x^ to A3. 1\) 
and that Assumption \4-l\ holds. Furthermore, assume that 

\\DF{x)\\ < L VxeB. (4.3) 

Let the positive stepsize, ji, be such that 

^ < q 



2G q Li 



1 G q Li .. ,\ ...^ 

2 q 



Let 



1+? 



(4.4) 



\(x ,x^)<p, (4.5) 
then the iterates satisfy 



A p (x k ,x^<p, A; = 1,2,... (4.6) 
and A p (x k ,x>) — > as k — > oo. Moreover, let 

c = C~^ (U-^lA. (4.7) 



,2 r q 
The convergence rate is given by 

A p {x k ,x^<p{l-c) k , (4.8) 

if e = 1. Fore e (0, 1), t/ie convergence rate is given by 

A p (x k ,x^) < (ckl^ + p- 1 ^) l£ , A; = 0,1,.... (4.9) 



Proof Using (12. 7p and (12.31) . we obtain, for the sequence of residues, 

= A p (x k ,x ] ) + - (||x fc+ i|| p - \\x k \\ p ) - (J p {x k+ i) - J p (x k ),x ] ) 

Q (4.10) 

= A p (x k ,x^) + - (|| J P (x k+1 )\\ q - \\J P (x k )\\ q ) 
-{J p {x k+ i) - J p (x fc ),x t ). 
Applying ( 12 .TP and (f) of Theorem 12.51 with x* = J p (xfc + i) and x* = J p (x k ), we get 

-(iij^+or-nj^or) 

g G (4.H) 

< — \\J p (x k+1 ) - J p {x k )\\ q + (J p (x k+1 ) - J p (x k ),x k ). 
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Substituting (12.121) and using this inequality in f )4.10p yields 

A p (x k+ i, art) - A p (x fc , art) 

< f\\vDF(x k )*j p (F(x k )-y)\\« (4.12) 
- (fiDF(x k )*j p (F(x k ) -y),x k -art). 

We estimate each term in ( 14. 12|) separately. The first term satisfies the estimate 

^DF{x k y 3p {F{x k )-yW < % g L%F(x k )-y)r. (4.13) 
q q 

For the second term, we have that 

-{/j l DF(x k )*j p (F(x k ) - y),x k - art) 
= -v(j P (F(x k ) - y), DF(x k ){x k - art)) 

= -K(jp( F ( x k) - y), F(x k ) - y) 

-(j p (F(x k ) - y), F(x k ) -y- DF(x k )(x k - art))). 

Note that, by the fundamental theorem of calculus for Frechet derivative, we obtain 
that 

\\F(x k )-y-DF(x k )(x k -x^)\\ < ^\\x k -x^\\ 2 . (4.15) 

Then, using (I2.8P and stability (c) of Assumption 14.11 we have 
-(nDF(x k )*j p (F(x k ) - y), x k - ar f ) 

= -fi\\F(x k ) - y\\ p 

+v{j P (F(x k ) -y),F(x k ) -y-DF(x k )(x k -art)) 

< -fi\\F{x k ) - y\\ p + ^L\\(F(x k ) - y)\r l \\x k - x^f 

( \ 2//p 

< ^iijr^-yllP+lLca^j ||F(x fc )-y|r £ . 

Combining these estimates and using the notation 

7 fc = A p (x k ,x j ), 

we obtain 

7 fc+ i-7.< (% 9 L'-±v) \\F{x k )-y\\> 

- -v\\F(x k ) - y\\ p + ^LCl [-^J \\F(x k )- y \\^. 

We claim that 

7 fc+ i = Ap(x fc+ i,ar t ) < p, (4.18) 
which we prove by induction. Assume that 

A p (ar m ,xt)<p (4.19) 
holds for m = 0,1, . . . ,k. With the mean value inequality, it follows that 

ll*W> - VW < I' ' = L(? MC P )^ ™ = 0, 1, 2, ■ ■ ■ , *■ (4-20) 
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Therefore, 

- ifi\\F(x m ) - y\\ p + ^LC 2 F (p/C p ) 2/p \\F(x m ) - yf +e < 0, (4.21) 
m = 0, 1, 2, . . . , k. Dropping this non-positive term, we obtain 

7fc+i - Ik < (y/^ 9 ~ \t^j \\F(x k ) - y\\ p . (4.22) 

Note that the term ( — -[i q L q — — juj ||F(xfc) — is non-positive. We obtain that 

A p {x k+1 ,x*) <p, (4.23) 

which establishes the claim. 

Now, we return to (I4.22p . By the Holder type stability (14. 2p . we have that 

lk+i <lk~ cryf* (4.24) 

Note that, by the conditions on fi, we have < c. By letting k go to infinity on both 
sides of the above inequality, we conclude that 

7fc — > as k — > oo. 

In the remainder of the proof, we obtain the convergence rate. Note that, with the 
choice ( 14. 4 ft of /i, 

< c < 1. (4.25) 

With e = 1, we have 

lk+i < (1 - c) 7fc (4.26) 

which expresses the convergence rate (14. 8ft . 

For the convergence rate with e G (0, 1), from ( I4.24p . we obtain that 

7j#>7^(l-C7^)-^. 



Noting that 



'l-xH+i>l + - — -x VxG(0,l), 
~ 1 + e 



we have that 
It follows that 



l=£ _i^£ 1-e 

l+s \ 1+e ■ x c 



1+e . . 1+e 

1 - £ ^ /1-e i- E 



For the critical index e = 0, we obtain 
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Theorem 4.6 Let Y be a general Banach space, and X be a Banach space which is 
p-convex and q-smooth with l/p + 1/q — 1. Assume there exists a solution x^ to A3. 1\) 



and that Assumption \4.S\ holds with e = 0. Furthermore, assume that 

\\DF(x)\\<L VxG#, (4.27) 
and that the stability constant Cp and the positive stepsize \i satisfy the inequality 



G q Li \ 2 * \C P 
Then the iterates satisfy 



2 - 



^<^- \1-\lC 2 AM 1. (4.28) 



A p (xfc,x"f) — )■ as k — >■ oo. 



Moreover, let 



c = fi[ % q -^ - 1 + \LCl [£\ ' ) c£*. (4.29) 



q 2 * \C P/ 

The convergence rate is given by 

A p (x fe , x f ) < (A p (x , xt)- 1 + cA;)- 1 . (4.30) 



Proof Using ( I4.17P in the proof of Theorem 14.51 subject to the substitution e = 0, we 
obtain that 



Tfe+i - Ik < f^V^ 9 ^LC F (jQ ||F(x fe ) - yf. 



Note that, by (I4.28p . the right-hand side of the above inequality is non-positive and 
< c < 1. Then, using the Holder type stability ( 14 .2p with e = 0, we have that 

lk+i <lk~ or*- ( 4 - 31 ) 

The convergence result and convergence rate (14.301) can be deduced by using the same 
arguments as in the proof of Theorem 14.51 



Remark 4.7 The Holder type stability condition A3.3\) or A,4-^ implied by a lower 
bound of the Frechet derivative DF . More precisely, if, there exists a constant C such 
that 

\\DF{x) ( u X ~ X l u ) || > C\\x - x+H 1 -" Vx G BJx*) n V(F), 
\ \\x — X'\\ J 

for some a G (0, 1] and r sufficiently small, then, by combining this and 

\\F(x) -F(x) - DF(x)(x - x)\\ < ^||x-xf Vx, x G V(F) 

we obtain that 

\\x - x f || < C F \\F{x) - F(x ] )\\^ Vx G B r (x^) n V(F), 
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for some constant Cp depending on C and L. The ill-posedness of many inverse problems 
indicates that in general it is impossible to obtain a lower bound for DF. If one 
projects the forward operator F properly, an estimate for the lower bound of DF could 
be obtained. Under various conditions, the lower bound for DF has been investigated in 
the analysis of inverse problems. For example, see Calderon [13], Somersalo, Cheney & 
Isaacson JEM/, and Dobson J7?| / for the electrical impedance tomography problem, Bao, 
Chen & Ma J2J/ for the inverse medium problem associated with the Helmholtz equation, 
and Ammari & Bao ^ for the inverse medium problem for electromagnetic waves. 

5. Example: Electrical Impedance Tomography 

In this section, we discuss Calderon's inverse problem, which forms the mathematical 
foundation of the Electrical Impedance Tomography (EIT) problem [13]. For a recent 
review, we refer to Uhlmann [34J. We mention some key uniqueness results, namely, 
by Kohn & Vogelius [231 121], Sylvester & Uhlmann [33], and Astala & Paivarinta [8]. 
Here, we focus on results pertaining to stability; see Alessandrini [3], HJ [5] . In particular, 
we relate to the work of Alessandrini & Vesella [6] and Beretta & Francini [10], who 
establish a Lipschitz type stability estimate if the conductivity is piecewise constant 
on a finite number of subdomains with jumps, for the real-valued and complex-valued 
cases, respectively. 

5.1. The Dirichlet-to-Neumann map 

Let Q C M n be a bounded domain with smooth boundary. The electrical conductivity 
of Q is represented by a bounded and positive function j(x). Given a potential 
/ G H l l 2 (dVL) on the boundary, the induced potential u G if x (f2) solves the Dirichlet 
problem 




0= o, 

«= f, 



in f2 
on d£l. 



The Dirichlet-to-Neumann map, or voltage-to-current map, is given by 




where v denotes the unit outer normal vector to dQ. 
The forward operator F is defined by 

F : Ici™(H)4 C{HV 2 {dty,H-V 2 {dn)), 

The Frechet derivative DF of F at 7 = 70 is given by 

DFfro) : IcL°°(n) C{Hy 2 (dn),H- 1 / 2 {d^l)) 
8l ^DF{^){5i), 



(5.1) 



(5.2) 
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and DF(j )(Sj) is given by 

(DF( l0 )(5j) f,g)= [ 5 7 V M • V^dx, f,ge H x ' 2 {dn) (5.3) 
Jn 

where 

V ■ (70 Vu) = V • (7 Vw) = 0, in ft, 
u — f, v = g on dQ. 

We note that L°°(n) is not a uniformly convex Banach space. Furthermore, to get the 
Holder type stability, the preimage space needs to be reduced. We specify the proper 
space X in Subsection 15.31 

For n — 2, Astala and Paivarinta proved that A 7 uniquely determines 7 under the 
assumption that 7 G L°°(Q). For n > 3, Paivarinta, Panchenko and Uhlmann [28] 
proved the uniqueness under the assumption that 7 G W 3 ^ 2,00 (fl). 

5.2. Lipschitz stability 

It is possible to obtain Lipschitz type stability, essentially, by assuming that 7 belongs 
to a particular finite dimensional space. 

We write x = (x',x n ) where x' G M™ -1 for n > 2. With B R (x), B R (x') and Qr(x) 
we denote respectively the open ball in K n centered at x of radius R, the ball in R™ -1 
centered at x' of radius R and the cylinder B' R {x') x [x n — R, x n + R). For simplicity of 
notation, -Br(O), B' r (0) and Qr(0) are denoted by Br, B' r and Qr. 

Definition 5.1 Let Q be a bounded domain in 1R™. We say that dQ is of Lipschitz 
class with constants ro, L > if, for any P G dQ, there exists a rigid transformation of 
coordinates such that P = and 

Q n Q ro = {{x',x n ) G Q ro I x n > 4>(x')} 

where <fi is a Lipschitz continuous function on B' ro with 0(0) = and 

U\\c°-i(B> ro ) < Lr . 

Definition 5.2 Let Q be a bounded domain in ~R n . Given a G (0,1), we say that 
dfl is of C 1,a class with constants r , L > if, for any P G dQ, there exists a rigid 
transformation of coordinates such that P = and 

Q n Q ro = {(x',X n ) G Q ro I X n > 4>(x')} 

where is a C 1,a function on B' ro with 0(0) = |V0(O)| = and 
UWc^(B') < Lr . 
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Assumption 5.3 Q C M. n is a bounded domain satisfying 
\Q\ < A\B ro \. 

Here and in the sequel |Oj denotes the Lebesgue measure ofQ. We assume that dQ is 
of Lipschitz class with constants r and L. 

Assumption 5.4 The conductivity 7 is a piecewise constant function of the form 

N 
3=1 

satisfying the ellipticity condition 
K~ l < 7 < K 

for some constant K , where = 1, . . . N are unknown real numbers and Dj are known 
open sets in M. n . 

Assumption 5.5 The Dj,j = 1,...,N are connected and pairwise non- overlapping 
open sets such that U^ =1 Dj = Q and dDj are of C 1,a class with constants ro and L 
for all j = 1,...,N. We also assume that there exists one region, say D\, such that 
dDi fl dQ contains an open portion Si of C 1,a class with constants ro and L. For every 
j E {2, . . . , iV} there exist ji, . . . ,jM £ {!)■■■) N} such that 

D J1 =D 1 , Dj M — Dj 

and, for every k = 1, . . . , M, 

dD jk _, n dD jk 

contains a open portion of C l,a class with constants ro and L. 

Alessandrini and Vessella [6] establish the following Lipschitz stability estimate 

Theorem 5.6 Let satisfy Assumption 15.31 and ^ k \k = 1,2 be two real piecewise 
constant functions satisfying Assumption \5.4\ and Dj,j = 1,...,N satisfying 
Assumption \5.5[ Then there exists a constant C = C(n, ro, L, A, , K, N) such that 

|| 7 (!) _ 7 ( 2 )|| Loo(n) < C||A 7 (1) - A 7 (2) ||£(tfl/2( 9 n),#-l/2(£>fi))- (5.4) 
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5. 3. Convergence 

We verify that the assumptions of Section H] can be satisfied. We specify our preimage 
space as 



equipped with L p -norm, p > 1. With the aid of this particular basis of X, one can 
show that F and DF are Lipschitz continuous. Moreover, assuming that 71,72 satisfy 
Assumption 15.41 and Q satisfies Assumption 15.31 we have the following estimates: 



where C, L and L depend on f2, N and ellipticity constant K. Furthermore, since X is 
finite dimensional, the weak topology is equivalent to the strong topology. Hence, F is 
a weakly sequentially closed operator. 

Let Q satisfy Assumption 15.31 preimage space X be defined by ( 15. 5ft and F be 
defined by ( 15. ip . Assume that y = F(j*) for some 7^ G X. Then Assumption 14.11 and 
( 14. 3 p of Theorem 14.51 are satisfied. Hence the Landweber iteration f |2.12j) converges with 
convergence radius given by ( 14.51) and convergence rate given by ( 14. 8h . Convergence 
of a regularized Newton method for a finite dimensional EIT problem was proven by 
Lechleitner and Rieder [26]. Their analysis, however, is based on the tangential cone 
condition. 

6. Discussion 

We discuss a Landweber iteration method for solving nonlinear operator equations in 
both Hilbert and Banach spaces. Traditionally, the gradient-type methods are often 
regarded as too slow for practical applications. Provided that the nonlinearity of the 
forward operator obeys a Holder type stability, we could prove the convergence and 
give a sublinear convergence rate. With a Lipschitz type stability, the convergence 
rate switches to a linear one. Based on these convergence rates, we anticipate that 
this Landweber iteration is a valuable tool in solving inverse problems in both Hilbert 
and Banach spaces. This also motivates the study of Holder /Lipschitz type stability in 
inverse problems to provide explicit reconstructions. 
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X = span{xD i; • • • ,Xd n } 



(5.5) 



||F(7x) - F(-f 2 )\\ C (m/2(n),H-m(n)) < c hi ~ 72|U*(n), 

\\DF\\ jC ^ x ,C(H 1 / 2 (Q),H- 1 / 2 (a))) < L, 

DFfa) - DF(-f 2 )\\ c{Hl/ 2 (nhH -i/2 m < LWjx - Tallin), 



(5.6) 
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